Dark states in the magnetotransport through triple quantum dots 
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We consider the transport tlirougli a system of tliree coupled quantum dots in a perpendicular 
magnetic field. At zero field, destructive interference can trap an electron in a dark state — a 
coherent superposition of dot states that completely blocks current flow. The magnetic field can 
disrupt this interference giving rise to oscillations in the current and its higher-order statistics as 
the field is increased. These oscillations have a period of either the flux-quantum or half the flux- 
quantum, depending on the dot geometry. We give results for the stationary current and for the 
shotnoise and skewness at zero and flnite frequency. 
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The quantum-mechanical interference of electronic 
paths in a conductor gives rise to a number of inter- 
esting phenomena in mesoscopic physics. Perhaps the 
most familiar is the occurrence of Aharonov-Bohm (AB) 
oscillations [H in the current through multiply-connected 
structures in a magnetic field These oscillations 

arise due to the accumulation of a phase difference 
between different paths through the device given by 
(f) ^ § A ■ d\ = 27r$/$o with (f> the flux enclosed by 
the device and = h/e the magnetic flux quantum. 

The period of such oscillations is dependent on the 
nature of the interfering paths, and therefore on the spe- 
cific system in question. Flux periods of $0 are what one 
expects in conventional AB experiments, such as those 
on coherent beams of electrons in free space This 
period is also frequently encountered in mesoscopic ex- 
periments, for example, in normal metal rings M and in 
electronic Mach-Zehnder-style interferometers [sl, Q, in- 
cluding those with one 0,11] or two Q quantum dots 
in the arms. The flux-period <&o/2 is also observed, not 
just in superconducting systems [lo|. but also in normal 
metals [ll|, [T^ due to weak- localisation effects [l^ . 

A different quantum coherent effect was described for 
mesoscopic systems in Ref. [3 — that of coherent pop- 
ulation trapping (CPT) in quantum dots. In this all- 
electronic analogue of a quantum optics effect 
the coupling geometry of a triple quantum dot (QD) leads 
to the establishment of a so-called "dark-state" that com- 
pletely blocks the current through the device. The dark- 
state is composed of a coherent superposition of elec- 
tronic states in different dots. 

In this paper we consider the interplay of coherent pop- 
ulation trapping and the AB phase. We demonstrate how 
a magnetic field can destroy the delicate phase-cancelling 
that maintains the dark-state, lift the current blockade 
and give rise to current oscillations as the field increases. 
Furthermore, we show that in a triple-QD structure such 
as in FiglTJ the oscillations can exhibit periods of both $0 
and ^$0 depending on the symmetry of the system. We 
give results not only for the stationary current but also 
for the shotnoise and skewness (second and third current 
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FIG. 1: Three quantum dots are coupled coherently to one 
another via tunnel couplings tij, and incoherently to source 
and drain leads with rates Fi. Each dot contains a single 
level and by adjusting the relative positions of these levels, 
the system can be prepared in a dark state where no current 
flows despite the applied bias. In a perpendicular magnetic 
field, the structure encloses a magnetic flux <I> which causes 
a phase-difference between different paths around the system 
that can disrupt the dark state and lead to current flow. 



cumulants, respectively) both at zero [iSj and at finite 
frequency • We also consider the effect of dephasing 
on the current oscillations. 

We consider the QD geometry depicted in Fig. [TJ The 
AB effect requires that all three dots be coupled in a ring 
structure as shown, which is in contrast to Ref. 



IJ] where 



only two 'bonds' were present, and no such effect would 
be observed. We work in the strong Coulomb blockade 
regime such that there is at most one excess electron in 
the three dot system at any one time. Each dot has a 
single level relevant to transport, and we denote as \i) 
the state with an electron in dot i . The Hamiltonian in 
the basis {|1), |2), |3)} is then 
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where e and A describe the energies of the dot levels and 
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tij are the tunnel couplings. Without magnetic field, 
time-reversal symmetry means that all Uj are real, and 
we take them all to be positive. Application of the mag- 
netic field breaks this symmetry and the amplitudes will 
be complex in general. We choose a gauge such that the 
phase (j) = 27r$ / $o is accumulated on the bond between 
dots 1 and 2. Finally, dots 1 and 2 are connected to 
source leads, and dot 3 to the drain. 

The density matrix (DM) for the system p{t) contains 
entries not only for the three single electron states 
but also the empty state |0). Within the Born-Markov 
and infinite-bias approximations, the time-evolution of 
the DM is given by the generalised master equation in 
the Lindblad form, 



dp 



= -t[H, p] + Y, D.pDl - \dId,p yolD^, (2) 
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where the quantum jump operators Di — v'T\|l)(0|, 
D2 = Vr^|2)(0|, and = Vr^|0)(3| describe irre- 
versible tunnelling of electrons into and out of the system 
with rates F^. In the following we set all these rates to 
equal: F^ = F. 

Starting from initial DM p(0) = |0)(0|, the DM 
at subsequent time has ten non-zero elements if we 
assume the most general parameters in Ti. These 
elements we arrange into the column vector pit) = 

(POO, Pll, P22, P33, 5P12, 5P13, 9P23, 5Jpi2, 3?Pl3, 5ip23) ■ 

The master equation ([2]) can then be written 



p = Mp, 



(3) 



with the time-evolution matrix M given in the Appendix. 
The stationary properties of the system are determined 
by the eigenvalues and eigenvectors of m[i^. The sta- 
tionary DM, p{oo) is given by the eigenvector of M with 
zero eigenvalue, whence the stationary average current 
(/) = Fp33(oo). To calculate higher-order statistics we 
require the full spectral decomposition of M [l^ . We ini- 
tially consider the behaviour of the system in the absence 
of dephasing and return its effects subsequently. 



ZERO-FIELD CURRENT 

We first consider the properties of the system at zero 
= 0). 



magnetic field 



In Ref. 



14l |. it was shown that, 
with the special choice of parameters ^12 = 0, e = A = 
and = t23, the system always reaches the trapped 
state l^*) = (|1) — |2)) /^/2 in the stationary limit (in 
the absence of dephasing). This exact parameter set is 
unlikely to pertain in experiment, and it is important 
to show the existence of the dark state for more general 
parameters. 

If we set A = Aq with 
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FIG. 2: At zero-field, the stationary current (7) through 
the three-dot system shows a pronounced anti-resonance with 
complete current blocking at a detuning of A = Ao where the 
dark-state forms. The three curves show the current for differ- 
ent values of the coupling parameters (ti2,t23), with tis = F 
in each case. At zero-field, a dark state always exists when- 
ever both ti3 and t23 are finite. 



then, as can easily be verified, the state 
1 



dark/ 



(i23|l)-<13|2)) 
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is an eigenstate of Hamiltonian ([T]). Moreover, the vector 
corresponding to the pure DM 



Pdark — | ^dark) (^dark | 



(6) 



is the eigenvector of M with eigenvalue zero, and thus 
Pdark IS the Stationary state of the system. Since this 
state has no electronic density at the drain dot (dot 3), 
the stationary current through the device is exactly zero. 
If we assume therefore, that we have experimental control 
of the detuning A by, for example, backgates under the 
dots, then a dark-state always can always be found at 
zero field by sweeping A. It should be noted that the 
detuning e does not affect the existence of the dark-state. 

A special instance of our geometry is when ti^ = t23 = 
to. In this case, the dark-state forms at a value of the 
detuning A = Aq = 0. With this choice, the system is 
symmetric under the exchange of dots 1 and 2, and cor- 
respondingly we will refer to this situation as the 'sym- 
metric' case in what follows. This special case admits a 
number of exact solutions. 

In this symmetric case at zero field with e = and A 
a freely variable parameter, we have for the stationary 
current 
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4{tl-tl,)^ + A^r^ + 6tl + 8tl,)- 



Figure [2] shows the stationary current through the device 
for B = as a function of the detuning A. We show not 
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only this result for the symmetric case, but also numerical 
results for the current with various different couplings. In 
each case an anti-resonance occurs with complete current 
suppression at A = Aq. In the following we will always 
set A = Aq, unless otherwise stated, such that there is 
CPT at zero-field. 

CURRENT OSCILLATIONS 

The application of magnetic field has the capacity to 
lift the dark-state current blockade and give rise to oscil- 
lations in the current. In the symmetric case, an exact 
expression for the current at finite field can be found. 
With e = 0, we have 

4(tf2+0+*?2(r2-*^(l + 7cos2</)))' 

which is plotted in Fig. [S] The current shows clear oscil- 
lations with superconducting flux period 5$o- We will 
return to a discussion of which paths interfere here in a 
moment. Let us first note that the maximum current as 
a function of $ occurs at $ = j^o, irrespective of ti2, 
and that the current at this value of $ is itself maximised 
by setting ti2 = to. With this choice, the current in the 
weak coupling limit to ^ T is given by (/) = 'ito/T sin^ 0, 
and in the opposite regime, ^ T, we have (/) = 2T/7 
for all (f) except at mr where it is exactly zero. This 
means that the oscillations are easier to observe in the 
weak coupling limit. 

The suppression of the current at $ = ±^$0, with n 
odd, is a consequence of the additional symmetry of the 
above situation. Figure [3] also shows the behaviour of 
the current as we move away from the symmetric cou- 
pling. As the asymmetry increases, the features at odd 
multiples of ^$0 disappear, doubling the period of the os- 
cillations to $0- Setting ti2 = t23 and ti^ = <23 (1 + a), 
a perturbation series for the current in a away from the 
symmetric case, shows the current at $ = i(f>o to be 

(/)-164aVr + 0(«4). (8) 

That this dependence is quadratic suggests that some 
degree of current suppression at odd- multiples of i$o 
may be visible in experiment. 

The above results can be understood by considera- 
tion of the interference between different paths around 
the dots. Consider the three dot system to be occupied 
and isolated from the leads. If we assume that at time 
T = the system is in the pure state |4'(0)), then the 
wavefunction at later time r is |*(t)) = e"''^*|vI/(0)) 
which, for small times, can be expanded as |^'(t)) ~ 
(1 - irH - t'^H'^) |«'(0)). Consider the system initially 
in the state |1). To first order is t, evolution under the 
full Hamiltonian with A = Aq gives 7i|l) = Ao|l) -I- 
ti2e~"^|2)-|-ti3|3). Therefore the first-order amplitude for 
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FIG. 3: The stationary current through the three-dot system 
shows pronounced oscillations as a function of the applied 
flux. Plotted here is 7 = (/)/{/)max, the ratio of the current 
to its maximum value, which always occurs at "J> = ^^o- 
The different curves are for different values of tia. Other 
parameters were ti2 = t23 — T, e = and A — Aq, such that 
CPT occurs at zero field. In the symmetric case with tis = t23 
(solid curve), CPT trapping occurs at f $o, n = 0, 1, 2, . . . and 
the flux-period of the oscillations is thus ^^o. As trA moves 
away from symmetry, the dark-state current blocking at ^$0, 
n odd disappears, and the period of the oscillations doubles 
to <I>o. 



the transmission from dot 1 to 3 is ag^"* = (3|7i|l) = ti^. 
Similarly, the amplitude from dot 2 to 3 is ~ ^23- 
Thus, if we start the system in the dark-state superpo- 
sition 1/N{t23\l) - ti3\2)), with norm N = \Jt\l+t\-i, 
these two paths interfere destructively at dot 3 with to- 
tal amplitude a^^^ = (^230^^ - ^I'iO^^hl^ = 0- The dark 
state is therefore stabilised against first-order tunnelling 
regardless of apphed field. 

The unblocking of the system at finite S-field occurs 
at second-order. Consider the second-order amplitude 

from 1 to 3: = (Aoii3 eti3 + e-''^ii2t23)/A^^ which 
has contributions from the three paths 113, 133, and 
123. Similarly, the three paths 223, 233 and 213 give 
the amplitude from dot 2 to 3 as = (— Aot23 + 
£^23 + e"^fi2ti3)/A^^. The total second-order amplitude 
a^"^^ for the dark-state electron to tunnel to dot 3 is 
a(2) ^ tx2 (e""^ - 1) (e"^t?3 -f /N'^- The correspond- 
ing probability = ja'^'p is 

= 2t\2 (1 ~ COS0) (4 + 4 + 2i?34 COS0) /N\ (9) 

We see immediately that this probability is zero at zero 
field. Furthermore, if ^ ^231 ^^e probability p^'^\ and 
hence the current, has period $o- However, if = t^^, 
this probability becomes p^^^ = 4tf 2^13 sin^ (j), which gives 
rise to current blocking with period 5<I>o. This halved 
oscillation period results from the symmetry of the two 
paths 123 and 213 when ^13 = t23. 
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This symmetry is reflected in the detuning required to 
obtain the dark state. When ^13 = ^23, we have Aq = 
= — Aq. However, when we have a couphng asymmetry, 
^0 7^ 7^ — Ag, and the period is doubled. Then setting 
A = — Ao, instead of A = Aq, is equivalent to shifting 
the phase of the current oscillations by 0, such that the 
dark-states occur at §$0; n odd. On the other hand, 
creation of a dark state at other values of <f> is impossible 
since this would require complex values of A, which is 
obviously a real parameter. 



DEPHASING 

We model the influence of dcphasing due to charge 
noise through the introduction of the three jump opera- 
tors 



D7 



* = 1,2,3, 



(10) 



which enter in Eq. ([2|) in the same way as do the jump 
operators Di above. Out of simplicity, we assume the 
decohercncc rate 7 to be the same for each dot and we 
give here only results for the case when all three couplings 
are equal, Uj = t. In this case, the current as a function 
of the phase (f> and decohercncc rate 7 is 



r(47(r + 27)t2 + 2/(0)) 



r7(r + 27)2 + 2(2r2 + i3r7 + 1472)^2 + 7/(^) 

with (/)-dependence contained in the function 



/(0) = 8i4(r + 37)(r + 27)- 



1 sin^ , 



(11) 



Since /(O) = ,the dephasing leads to a finite current 
at zero-field through the disruption of the coherence be- 
tween the two dark-state dots 141. The effect of de- 



phasing on the oscillations can be quantified through the 
visibility: 



1.(7) = 1 - 



(/)(0max,7) 



This is found to be 



4r(r + 3^)tl 



rj{r + 27)2 + 2t2(2r2 + i3r7 + 1472 



(12) 



(13) 



For small dephasing, 7 ^ F, the visibility deviates from 
unity as 



, 7 F 

"=^-^l2f + 4l 



(14) 



and for strong dephasing 7 ^ F, wc have ly = 3^^/72 . 




FIG. 4: The zero-frequency Fano factor F'^'(O) for tlie three 
dot system as a function of magnetic flux also exhibits os- 
cillations. We observe strong superPoissonian peaks with 
F''^\0) = 3 at the values of flux for which the dark-state 
forms [l^- Away from these points, the noise is subPoisso- 
nian. The period of these oscillations is the same as for the 
current. Same parameters as for Fig. [3] except for the dis- 
played values of ti3. 



HIGHER-ORDER CURRENT STATISTICS 

The dark state and magnetic field also leave their mark 
on the higher statistics of the current through the device. 
We consider here the two finite-frequency Fano-factors: 
F(2)(w) = S'(2)(a,)/(/) where S^^\uj) is the shotnoise, 
and {lu,lu') = S^^'' {uj,uj')/ {I) with the skewness 



5(3) (, 



^' {5I{0)SI{t)SI{t')).{15) 



These quantities can be straightforwardly calculated [l9j 
from Ai, the eigenvalues of Af, and V the corresponding 
matrix of its eigenvectors. Let Cj be the jump operator 
that transfers an electron to the drain from dot 3. In the 
basis of the vector p, it has elements {Cj)ij = TSuSji. 
The Fano factors can then be expressed solely in terms of 
the eigenvectors and the quantities Ck = {V~^CjV)kk- 
The finite-frequency shotnoise Fano factor is given by 



F(2)(^) = I-2Y: 



CkXk 



(16) 



The zero-frequency result, i^(2)(o), is shown in Fig. [4] 
as a function of for several couplings. At <!> = n^O: 
71 = 0, 1, . . ., we see highly superPoissonian maxima with 
Fano factor F(2)(0) = 3, which is the same value as 
found for the ii2 = model discussed in Ref. [l3|- In 
the symmetric case fi3 = t23, we see further maxima, at 
<i> = f $0, also with F(2)(0) = 3. These latter disappear 
as the coupling asymmetry increases. In between these 
sharp superPoissonian peaks, the shotnoise is strongly 
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FIG. 5: (colour online) Contour plot of the finite-frequency 
Fano factor F^'^\lj) as a function of magnetic flux <1> and fre- 
quency uj. Colours white, red, and blue correspond to Pois- 
sonian, super-, and sub-Poissonian values respectively. Large 
superPoissonian peaks occur in the zero frequency limit only. 
The shotnoise also shows a series of inflexion points at a set 
of frequencies corresponding to the energy differences AEij 
of the Hamiltonian (dotted lines). Same parameters as Fig. |3] 
with all tij — T . 



subPoissonian, further illustrating the dramatic variation 
in systems behaviour as the field is changed. 

The shotnoise Fano factor at finite frequency F'-^'>{uj) 
is shown for symmetric coupling in Fig. [5l We see that 
the large superPoissonian peaks occur only close to zero- 
frequency, as the majority of the behaviour is subPois- 
sonian. Nevertheless, further structure is to be observed 
at finite frequency, with a number of inflexion points oc- 
curring as a function of uj, the locations of which are 
determined by the spectrum of the isolated Hamilto- 
nian Ti. The inflexions are located at w = AEij, where 
AEij = \Ei ~ Ej\ are the differences between all the 
eigenenergies of H. In the case where all three couplings 
are equal, these energies are obtained from the three so- 
lutions of the equation Ef — St^Ei — 2ig cos (p ~ 0, and the 
corresponding differences are shown overlaid on Fig. O 

The skewness is calculated from an expression similar 
to Eq. ((T6|) . but lengthier. We obtain the zero- frequency 
result shown in Fig. [6l which has the zero-field limit of 
F(3)(0) = 13, which is again the same as in Ref . [3| ■ The 
behaviour of the skewness as a function of field strength 
is similar to that of the shotnoise, but here the contrast 
between the values with and without the dark start is 
even more pronounced. 

In Fig. [7] we plot the finite- frequency skewness for sev- 
eral values of magnetic field with symmetric couplings. 
For $ ~ f ^o; n = 0, ±1, ... in this symmetric case, the 
skewness shows a sharp superPoissonian peak at the ori- 
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FIG. 6: The zero-frequency skewness Fano-factor ^'"^^(0) as 
a function of flux The behaviour is similar to that of the 
shotnoise, with the highly superPoissonian maximum value of 
-F'^'(0)max = 13. Parameters as in Fig.H 



gin and also strong superPoissonian behaviour along the 
symmetry lines of F^^^^ . Away from these values of the 
magnetic field, the skewness is strongly subPoissonian, 
except for slight SuperPoissonian shoulders at interme- 
diate frequencies before the Poissonian value of unity is 
recovered at high frequency The fine structure in 

the skewness arises from resonances between the three 
frequencies |w|, |w'| and |a; — w'l and the energy differ- 



ences AEi. 



CONCLUSION 

The behaviour of the coupled triple quantum dot sys- 
tem in a perpendicular magnetic field studied here is gov- 
erned by the interplay of two quantum-coherent effects: 
coherent population trapping and the Aharonov-Bohm 
phase. 

We have shown that a dark-state, for which no current 
flows, exists at zero field for arbitrary couplings between 
the dots provided ti^ and i23 are both finite. The mag- 
netic field can lift destructive interference maintaining 
the dark state, and give rise to oscillations in the current. 
For arbitrary parameters the period of these oscillations 
is $0j but in the special case when the coupling strengths 
ti3 and t23 are equal the period is halved to ^^o- These 
oscillations are also visible in the zero-frequency shot- 
noise and skewness which show large oscillations between 
strong superPoissonian and subPoissonian behaviour. Fi- 
nally, at finite frequency these quantities show consider- 
able structure which again show dramtic dependence of 
the magnetic field. 
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FIG. 7: (colour online) Contour plots of the finite-frequency 
skewness Fano-factor F^^\lj) as a function of its two fre- 
quency arguments ui and ui' for values of the magnetic flux 
^/^o = 0, 1/12, 1/6, 1/4. Parameters and colour scheme as 

in FigO Strong superPoissonian behaviour only occurs close APPENDIX: TIME-EVOLUTION MATRIX 

to $/$o = f , otherwise the skewness is predominantly sub- 
Poissonian. The fine structure arises from resonance between 
the frequencies and |a; — tj'l and the level-splittings 

of the isolated dot-system. 

The equation-of-motion matrix for the system is 
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